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1. Introduction 

The supersymmetry (SUSY) is believed as a symmetry of the unification theory such as su- 
pers tring theory and supersymmetric gauge theory is a candidate of a theory beyond the Standard 
model. However, it is broken in our current universe anyway. Since it cannot be broken by higher 
loop effects in perturbation, it is important to study the breaking nonperturbatively. 

Witten index [jl|] is a useful index related to the spontaneous SUSY breaking, which is defined 
nonperturbatively. Using the fermion number operator F, it is given by the following trace, 

w = tr(-l) F e-? H = (N B -N F )\ E=Q , (1.1) 

where H is the Hamiltonian of the system and E is its eigenvalue. As long as the spectrum is 
discrete, the index does not depend on a parameter p. 1 It is simply a difference of numbers of 
bosonic supersymmetric vacua and fermionic vacua. If the index is not zero, there exists at least 
one supersymmetric vacuum so SUSY is not broken. But if the index is zero, SUSY may or may not 
be broken, since it can be a result of cancellation between bosonic and fermion vacua, or a result 
of no supersymmetric vacua at all. The purpose of this talk is to propose a method to measure 
the Witten index using lattice simulation based on Ref. [||]. For a different approach from lattice 
simulation, see Ref. [J3]] . 

In terms of the path integral, the index becomes a partition function with periodic boundary 
condition [@, §] 

w = Zp = y^0^v7^exp(-5 P ), (1.2) 

where <p is boson, y and \j7 are fermion, and subscript P stands for periodic boundary conditions 
for all the fields in the temporal direction. It seems difficult to measure this quantity using lattice 
simulation, since what we usually measure is an expectation value normalized by the partition 
function but we need the normalization factor here. The normalization of the path integral measure 
is relevant as well. 

In the following section, we will discuss how to obtain the correct normalization of the par- 
tition function and thus the Witten index. And then in section |3| we confirm that it in fact works 
in supersymmetric quantum mechanics of which the Witten index is well known using a lattice 
simulation. We also test a method which would improve the efficiency of the measuring the Witten 
index. 

2. Idea 

We have to determine two normalizations: one for the path integral measure and the other is 
for the partition function (from the lattice data). 

The path integral measure with a correct normalization is easy to obtain. We only have to 
follow a standard derivation of the path integral from the operator formalism, where we insert 
normalized complete sets at each of discretized time slices ( Fig. [I] ). Regarding the discretization 

'if the spectrum is continuum, one has to take a limit j8 — > °°. 
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T = aN 




insert complete sets 

Figure 1: Obtaining the correct measure: Derivation of the path integral is exactly the lattice regularization. 



is the lattice discretization, we obtain the following measures for bosons and fermions: 
Bosons: J 
Fermions: 
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(2.1) 
(2.2) 



The correct normalization of the partition function is non-trivial. Let us start with a 1 -dimensional 
bosonic system with N lattice sites and consider the following quantity: 



e +5 exp 



c 



(2.3) 



regularization functional 

where }i is an arbitrary (positive) real number which should be tuned later and 

1 



C 



^(j) exp 



;I> 2 (0i 
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-N 



(2.4) 



Here, we have used eq. ( |2. Combining eq. ( |2.3[ ) and ( |2.4[ ), we obtain 

C 



-s 



exp 



(2.5) 



Since we can calculate the value of C analytically, and the denominator in the r.h.s is an observable 
in the lattice simulation, we can measure the partition function Z. Notice that though we have used 



a gaussian functional as a regularization functional in eq. (2.3), one can use any functional as long 



as it gives a calculable and convergent value like in eq. (2.4). 



In the r.h.s. of eq. ([2.5|), the action S appears with a "wrong sign" which cancels the original 
distribution. That is, the partition function is calculated using an extreme reweighting. To obtain a 
better efficiency, we have to tune the value of /i. 

Next let us introduce fermions. After integrating out the fermions, we obtain the effective 
action as usual: 

5' = 5 B -lii|det(D)|, (2.6) 

where 5b is the bosonic part of the action and D is the the fermion bilinear operator (i.e., the 
Dirac operator plus the Yukawa interactions) 2 . The phase factor of det(D) should be reweighted 



"If the fermion is Majorana, the determinant should be replaced with a Pfaffian 
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afterwards which gives for arbitrary expectation values 

I^Aa[D]e- s ' _ (Aa[D]) Q 
W f®4>o[D]e-* (<J[D})o , { } 

where o[D] is the phase factor and the subscript stands for a phase quenched average. This 
time we have to cancel e a factor o[D]e~ s ' to obtain the partition function. Therefore, measuring 
(a[D] _1 e +5 'exp(-^,^ 2 (0/ at ) 2 )), we obtain the Witten index as 

(ct[»p])o,p 



<exp[+^-iL.At 2 0, 2 ] /Oi) 



w = Z P =C-, r , .\ , r r .^m > ( 2 - 8 ) 

where C is given in eq. ( |2.4| ). 



The r.h.s of eq. (2.7) implies the phase quenched average of the phase factor o[D\ is almost 



the partition function. This observation is correct, and eq. (2.8) provides the correct normalization 
to the partition function. 

3. Numerical Test: Supersymmetric Quantum Mechanics 

We test our method using supersymmetric quantum mechanics (of j¥ = 2 Wess-Zumino type) 
of which the Witten index is known. 

If the lattice action keeps a part of supersymmetry as an exact symmetry on the lattice, we 
expect that the Witten index is well defined. More precisely, if the action is given as QA with an 
exact supertransformation which satisfies Q 2 = 0, we can repeat a similar argument to the contin- 
uum case. As a result, the index is well defined even at finite lattice spacing in such lattice models. 
In particular, the index from a finite lattice spacing should be an integer. 

A Q-exact lattice action for the supersymmetric quantum mechenics is given as [0] 

5 = L " + \ w '^k? + - fe)W(fe) - \fu +W k (Vk + i - Wk) +W"(<!> k )y k y k \ , 

k=0 LZ L z J 

(3.1) 

where ^ is a real boson, \ffk an d Yk we fermions, and F% is a real bosonic auxiliary field. The 
potential W is a function of and the prime (') indicates a derivative. If the asymptotic behavior 
is oo) > the supersymmetry is not broken and W(+°°)W(— °°) < it is broken. We 

use the following two cases: 

• n = 4: W = A 4 4 + A 2 2 SUSY, w = \ 

• n = 3: W = + h<!> 2 ^SYT w = 

where A, are parameters of the potential. We use the Hybrid Monte Carlo algorithm. See [^] for 
the implementation for this system. 

The results are plotted in Figs. || and ||. With a suitable choice of the known indexes 
are reproduced. There is almost no dependence on the lattice spacing, as expected from the exact 
<2-symmetry of the action. 

at 

Next, we consider a possible way to improve the efficiency. Because of the factor e p in 



eq. (2.8), the efficiency is poor and we need large statistics. This factor cancels the weight from the 
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Figure 2: n = 4 case, where the index is known to be 1 . Lis the physical size of the system, (left panel) 
jl 2 dependence, (right panel) Lattice spacing a = 1 /N dependence, (bottom) values of [I 2 and the measured 
index w, which minimize the error. 
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Figure 3: n = 3 case, where the index is known to be 0. L is the physical size of the system, (left panel) 
jl 2 dependence, (right panel) Lattice spacing a = 1 /N dependence, (bottom) values of \i 2 and the measured 
index w, which minimize the error. 



action so we do not have to use importance sampling with respect to a weight factor e~ s . Therefore, 
we can also use configurations generated with less importance sampling. Decomposing the weight 
factor as e~ s = e~ rS e~( l ~ r ^ s , we rewrite a general expectation value as 

f^A g - r VC^ _ (Ae- rS ) r 
{ ' J^e-^e-i 1 -^ {e- rS ) r ' K ' 

where ( ■ ) r is an expectation value with a weight factor e 1 - 1- '"' 5 . Therefore, preparing configurations 
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using r ^ s ' p , we can obtain the Witten index as follow: 

{o{D P ]e- rS v} r>P 



w = C 



(3-3) 



Note that r = is the usual importance sampling. 

We plot the result from the less importance sampling in Fig. f|. On the left panel, we see 
that the correct index is reproduced with a suitable choice of /I 2 . On the right panel, we plot the 
behavior of the errors versus number of the configurations used in the measurements. Contrary to 
the naive expectation, the magnitudes of the error are the same for large statistics in both r = case 
and r > case. For small statistics, however, r > cases converge to a line (num. of confs.) -1 / 2 
faster than r = case. This implies that the less importance sampling method is robuster for small 
statistics. 
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Figure 4: Results from the less importance sampling, (left panel) The obtained Witten index. Set la- 
bels are the same in Figs. ^| and [}| (right panel) Behavior of the errors, for set 4b. The dotted line is 
(num. of conf.) -1 / 2 . 



4. Conclusion and Discussion 

We proposed a method for measuring the Witten index, which is a useful index to detect a 
spontaneous supersymmetry breaking. Since the index is given as a partition function under the 
periodic boundary condition, it is important to use the correct normalization of the path integral 
measure. We also normalized overall factor of the partition function measuring a special regular- 
ization functional. As a test of the method, we measured the index of supersymmetric quantum 
mechanics. The results reproduced the known values of the index. A disadvantage of the method 
is its poor efficiency. A less importance sampling method may improve it to some extend. 

Finally, we mention possible applications of the method, which may or may not be practical. 
It is straightforward to use the method in higher dimensional systems. Within one-dimensional 
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systems, the most interesting one is supersymmetric Yang-Mills quantum mechanics with 16 su- 
percharges. This model is one of the candidates of M(atrix)-theory, and assumes the Witten index 
should be 1 to obtain a suitable supergravity limit. 
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